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Abstract 

We show the asymptotic long-time equivalence of a generic power law 
waiting time distribution to the Mittag-Leffier waiting time distribu- 
tion, characteristic for a time fractional continuous time random walk. 
This asymptotic equivalence is effected by a combination of "rescal- 
ing" time and "respeeding" the relevant renewal process followed by 
a passage to a limit for which we need a suitable relation between 
the parameters of rescaling and respeeding. As far as we know such 
procedure has been first apphed in the 1960s by Gnedenko and Ko- 
valenko in their theory of "thinning" a renewal process. Turning our 
attention to spatially one-dimensional continuous time random walks 
with a generic power law jump distribution, "rescaling" space can be 
interpreted as a second kind of "respeeding" which then, again under 
a proper relation between the relevant parameters leads in the limit 
to the space-time fractional diffusion equation. Finally, we treat the 
'time fractional drift" process as a properly scaled limit of the counting 
number of a Mittag-Leffier renewal process. 
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1 Introduction 



The purpose of this paper is to outhne the fundamental role the Mittag- 
Leffler function in renewal processes that are relevant in the theories of 
anomalous diffusion. As a matter of fact the interest in this function 
in statistical physics and probability theory has recently increased as is 
shown by the large number of papers published since 1990 of which a brief 
(incomplete) bibliography includes [SliailSlEIlEaEIllMlEgEDlEIllMl 

El El in EH [601 EH EZ]. 

In this paper we develop a theory for long-time behaviour of a renewal 
process with a generic power law waiting distribution of order /?, < /? < 1 
(thereby for easy readability dispensing with decoration by a slowly varying 
function). To bring the distant future into near sight we change the unit of 
time from 1 to 1/r, < r ^ 1. 

For the random waiting times T this means replacing T by tT. Then, 
having very many events in a moderate span of time we compensate this 
compression by respeeding the whole process, actually slowing it down so 
that again we have a moderate number of events in a moderate span of 
time. We will relate the rescaling factor r and the respeeding factor a 
in such a way that in the limit r — > we have a reasonable process, 
namely one whose waiting time distribution is the Mittag-Leffler waiting 
time distribution whose density is 

cl>''''{t) = -j^Epi-t''), </?<!, (1.1) 
with the Mittag-Leffler function 

We will call the renewal process with waiting time density (j)^^{t) the 
Mittag-Leffler (renewal) process. This process can be seen as a fractional 
generalization of the Poisson process, see |40] . 

Our method is, in some sense, analogous to the one applied in the Sixties 
of the past century by Gnedenko and Kovalenko [TS] in their analysis of 
thinning (or rarefaction) of a renewal process. They found, under certain 
power law assumptions, in the infinite thinning limit, for the waiting time 
density the Laplace transform 1/(1 -|- s^) but did not identify it as a Mittag- 
Leffler type function. In Section 2, we provide, in our notation, an outline of 
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the thinning theory for renewal processes essentially following Gnedenko and 
Kovalenko. Their method has inspired us for the reatment of our problems. 

As we consider our renewal process formally as a continuous time random 
walk (CTRW) with constant non-random jumps 1 in space (for the counting 
function N{t), in Section 3 we embed ab initio our theory into that of the 
CTRW, thus being in the position to treat the theory of a time fractional 
CTRW as limiting case of a CTRW with power law waiting time distribution. 
In this context the pioneering paper by Balakrishnan [Ij of 1985 deserves to 
be mentioned. Balakrishnan already found the importance of the Laplace 
transform 1/(1 + s^) in the time fractional CTRW and time fractional 
diffusion, but also did not identify it as the Laplace transform of (f)^^^{t). 
Then, in 1995 Hilfer and Anton [29], see also [271 I28j . showed that this 
waiting time density is characteristic for the time fractional CTRW and can 
be expressed in terms of the Mittag-Leffler function in two parameters, that 
is 

^^'\t)=t''-^Ep^p{-tf), 0</3<l, (1.3) 
with the generalized Mittag-Leffler function 

oo „ 

Ep^,{z) :=X; r(/?I + 7) ' ^""^ ^^-^^ 

The form (1.3) is equivalent to the form (1.1) that we prefer as it exhibits 
visibly also the cumulative probability function, the survival function, 

We explain in Section 4 two manipulations, rescaling and respeeding and 
use these in Section 5 to deduce the asymptotic universality of the Mittag- 
Leffler waiting time density under a power law assumption for the original 
waiting time. Then, in Section 6, assuming a suitable power law also for the 
spatial jumps we show that by a rescaling of the jump widths by a positive 
factor h (that means a change of the unit of space from 1 to l//i to bring 
into near sight the far-away space) another respeeding is effected, now an 
acceleration, that in the limit h ^ (under a proper relation between h and 
r) leads to space-time fractional diffusion. 

In Section 7, we pass to a properly scaled limit for the counting function 
N{t) of a renewal process (again under power law assumption) and obtain 
the time fractional drift process (viewing N{t) as a spatial variable). 
We will extensively work with the transforms of Laplace and Fourier, so 
easing calculations and proofs of convergence in distribution (also called 
"weak convergence") for our passages to the limit. 
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Essentially, we treat in this paper three topics. First, in Section 2, the 
thinning of a pure renewal process. Second, in Sections 3-6, under power 
law assumption for the waiting time, the asymptotic relevance of the Mittag- 
Leffler law, and then the general CTRW with special attention to space and 
time transition limits to fractional diffusion. As a third topic, in Section 7, 
we investigate the long time behaviour of the Mittag-LefSer renewal process. 

Essential properties of the derivative of fractional order in time and in space 
are given in Appendix A and Appendix B, respectively. Finally, in Appendix 
C we give details on the two special functions of the Mittag-LefHer type that 
play a fundamental role in this paper, the Mittag-Leffler survival probability 
and the Mittag-Leffler waiting time density. 

2 An outline of the Gnedenko-Kovalenko theory 
of thinning 

The thinning theory for a renewal process has been considered in detail by 
Gnedenko and Kovalenko [15]. We must note that other authors, like Szantai 
[631 [64j speak of rarefaction in place of thinning. Let us sketch here the 
essentials of this theory: in the interest of transparency and easy readability 
we avoid the possible decoration of the relevant power law by multiplying 
it with a slowly varying function. As usual we call a (measurable) positive 
function a(y) slowly varying at zero if a{cy)/a{y) 1 with y ^ 0^ for 
every c > 0, slowly varying at infinity if a[cy)/a{y) — > 1 with y +oo for 
every c > 0. A standard example of a slowly varying function at zero and 
at infinity is | logy|'^, with 7 E R. 

Denoting by tn, n = 1,2,3,... the time instants of events of a renewal 
process, assuming = io < *i < ^2 < ^3 < • • with i.i.d. waiting times 
Ti = ti , Tk = tk — tk-i for k > 2, (generically denoted by T), thinning 
(or rarefaction) means that for each positive index k a decision is made: 
the event happening in the instant t^ is deleted with probability p or it is 
maintained with probability q = l— p, 0<q<l. This procedure produces 
a thinned or rarefied renewal process with fewer events (very few events if q 
is near zero, the case of particular interest) in a moderate span of time. 

To compensate for this loss we change the unit of time so that we still have 
not very few but still a moderate number of events in a moderate span of 
time. Such change of the unit of time is equivalent to rescaling the waiting 
time, multiplying it with a positive factor r so that we have waiting times 
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tTi,tT2,tTs, . . ., and instants rti, rt2, rts, . . ., in the rescaled process. Our 
intention is, vaguely speaking, to dispose on r in relation to the rarefaction 
parameter q in such a way that for q near zero in some sense the "average" 
number of events per unit of time remains unchanged. In an asymptotic 
sense we will make these considerations precise. 

Denoting by F{t) = P{T < t) the probability distribution function of 

the (original) waiting time T, by f(t) its density (/(t) is a generalized 
function generating a probability measure) so that F{t) = £f{t')dt', and 
analogously by Fk{t) and /fc(t) the distribution and density, respectively, of 
the sum of k waiting times, we have recursively 

hit) = fit) , fkit) = f fk-iit - t') dFit') , for > 2 . (2.1) 

Observing that after a maintained event the next one of the original process 
is kept with probability q but dropped in favour of the second-next with 
probability pq and, generally, n — 1 events are dropped in favour of the 
n-th-next with probability p^"^ q, we get for the waiting time density of the 
thinned process the formula 

oo 

gq{t)=Y,QP''''fnit). (2.2) 

n=l 

With the modified waiting time r T we have 

PirT <t)= P(r < t/r) = FH/t) , 

hence the density /(t/T)/r, and analogously for the density of the sum of 
n waiting times /„(t/r)/r. The density of the waiting time of the rescaled 
(and thinned) process now turns out as 

oo 

5<z,r(t) = E«^'""'/-(*/^)/^- (2-3) 

n=l 

In the Laplace domain we have /n(s) = ' hence (using p = 1 — q) 

from which by Laplace inversion we can, in principle, construct the waiting 
time density of the thinned process. By rescaling we get 

HAS) = f;,p"-' (/>,))" = f^'\ . (2.5) 
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Being interested in stronger and stronger thinning (infinite thinning) let us 

now consider a scale of processes with the parameters r (of rescaling) and 
q (of thinning), with q tending to zero under a scaling relation q = q{T) yet 
to he specified. 

We have essentially two cases for the waiting time distribution: its 
expectation value is finite or infinite. In the first case we put 

A = / t' f{t')dt' < oo. (2.6a) 
Jo 

In the second case we assume a queue of power law type (dispensing with a 
possible decoration by a function slowly varying at infinity) 

^{t):= f{t')dt' ^ -t-^ ,t^oo if 0</3<l, (2.66) 

Then, by the Karamata theory (see \i2\ 168] ) the above conditions mean in 
the Laplace domain 



f{s) = l-Xs'' + l^s^j , for s ^ 0+ , (2.7) 

with a positive coefficient A and < (3 < 1. The case (3 = 1 obviously 
corresponds to the situation with finite first moment (2.6a), whereas the case 
< /? < 1 is related to a power law queue with c = XT{j3 + 1) sin(/37r)/7r . 
Now, passing to the limit of g — > of infinite thinning under the scaling 
relation 

q = XT'^, 0</3<l, (2.8) 

between the positive parameters q and r, the Laplace transform of the 
rescaled density g^{s) in (2.5) of the thinned process tends for fixed s 
to 

which corresponds to the Mittag-Leffler density 

g{t) = -±Ef,{-t0) = cl^''^{t). (2.10) 

Let us remark that Gnedenko and Kovalenko obtained (2.9) as the Laplace 
transform of the limiting density but did not identify it as the Laplace 
transform of a Mittag-Leffler type function. Observe that in the special case 
A < oo we have [3 = 1, hence as the limiting process the Poisson process, as 
formerly shown in 1956 by Renyi 
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3 The continuous time random walk (CTRW) 



The name continuous time random walk (CTRW) became popular in physics 
after Montroll, Weiss and Scher (just to cite the pioneers) in the 1960's 
and 1970's pubHshed a celebrated series of papers on random walks for 
modelling diffusion processes on lattices, see e.g. [Ml US], and the book by 
Weiss [6,6j with references therein. CTRWs are rather good and general 
phenomenological models for diffusion, including processes of anomalous 
transport, that can be understood in the framework of the classical renewal 
theory, as stated e.g. in the booklet by Cox In fact a CTRW can be 
considered as a compound renewal process (a simple renewal process with 
reward) or a random walk subordinated to a simple renewal process. 

A spatially one-dimensional CTRW is generated by a sequence of in- 
dependent identically distributed (iid) positive random waiting times 
Ti, T2, T3, . . . , each having the same probability density function (/>(t) , t > , 
and a sequence of iid random jumps Xi, X2, X3, . . . , in R , each having the 
same probability density w{x) , x G R. 

Let us remark that, for ease of language, we use the word density also 
for generalized functions in the sense of Gel'fand & Shilov [13], that can 
be interpreted as probability measures. Usually the probability density 
functions are abbreviated by pdf. We recall that (/)(t) > with (j){t) dt = 
1 and w{x) > with w{x) dx = 1. 

Setting to = , t„ = Ti -|- T2 -|- . . . T„ for n G N , the wandering particle 
makes a jump of length X„ in instant t„, so that its position is xq = 
for < t < Ti = ti , and x„ = Xi + X2 + . ■ ■ Xn , for t„ < t < tn+i ■ 
We require the distribution of the waiting times and that of the jumps to 
be independent of each other. So, we have a compound renewal process (a 
renewal process with reward), compare [7]. 

By natural probabilistic arguments we arrive at the integral equation for the 
probability density p(x, t) (a density with respect to the variable x) of the 
particle being in point x at instant t , see e.g. [211 ESI HH EHl [59l [60]. 



p{x,t) = 5{x)^{t) + I <p{t-t') 
Jo 



w{x — x') p{x' , t') dx' 



dt' , (3.1) 



in which the survival function 



(pit') dt' 



(3.2) 
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denotes the probability that at instant t the particle is still sitting in 
its starting position x = 0. Clearly, (3.1) satisfies the initial condition 
p{x, O"*") = 6{x). 
Note that the special choice 

w{x) = 5{x - 1) (3.3) 

gives the pure renewal process, with position x{t) = N{t), denoting the 
counting function, and with jumps all of length 1 in positive direction 
happening at the renewal instants. 

For many purposes the integral equation (3.1) of CTRW can be easily treated 
by using the Laplace and Fourier transforms. Writing these as 

oo 



^ />oo 

C{fit);s} = fis):= / e-'*fit)dt, 
Jo 



hoo 



J^{gixy,^}=giK):= I e+"^'' gix) dx , 

J — oo 

then in the Laplace- Fourier domain Eq. (3.1) reads 

p(k,s) = ^^^^ + <P(s)w(k)p(k,s) , (3.4) 

s 

Introducing formally in the Laplace domain the auxiliary function 

^(s) = = iifl , hence ^{s) = , (3.5) 

s^{s) 0(s) l + sH{s) 

and assuming that its Laplace inverse H{t) exists, we get, following [H], in 
the Laplace-Fourier domain the equation 

H{s) sp{k, s) — 1 = [Q{k) — 1] p{k, s) , (3-6) 

and in the space-time domain the generalized Kolmogorov- Feller equation 

/ H{t — t') —p{x,t') dt' = —p{x,t) + I w{x — x)p{x ,t) dx' , (3.7) 
Jo Ot J_oo 

with p{x, 0) = 6{x). 

If the Laplace inverse H{t) of the formally introduced function H{s) does 
not exist, we can formally set K{s) = 1/H{s) and multiply (3.6) with K{s). 
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Then, if K{t) exists, we get in place of (3.7) the alternative form of the 
generalized Kolmogorov-Feller equation 

d /"* r f~^°° 1 

—p{x,t) = / K{t-t') -p{x,t')+ / w{x-x')p{x',t')dx' dt' , 

(3.7') 

with p(.x, 0) = 5{x). 

Special choices of the memory function H{t) are (i) and (ii), see eqs (3.8) 
and (3.12): 

(i) H{t) = 5{t) corresponding to H{s) = 1 , (3.8) 
giving the exponential waiting time with 



(3.9) 



In this case we obtain in the Fourier- Laplace domain 

sp{k,, s) — 1 = [w{k) — 1] p{k, s) , (3.10) 
and in the space-time domain the classical Kolmogorov-Feller equation 

—p{x,t) = —p{x,t)+ / w{x — x') p{x' ,t) dx' , p{x,0) = S{x) . (3.11) 

(^t J —oo 



(ii) Hit) 



giving the Mittag-Leffler waiting time with 



, < /? < 1 , corresponding to H{s) = s^ ^ , 

(3.12) 



d 



(3.13) 

In this case we obtain in the Fourier-Laplace domain 



sp{k, s) — 1 = [w{k) — 1] p{k, s) . 



(3.14) 



and in the space-time domain the time fractional Kolmogorov-Feller equation 

/ + 0O 
w{x- x')p{x',t)dx' , p{x,0+) = S{x) , 
-oo 

(3.15) 



where tD^ denotes the fractional derivative of of order /? in the Caputo 
sense, see Appendix A. 

The time fractional Kolmogorov-Feller equation can be also expressed via 
the Riemann-Liouville fractional derivative tD^~^ , see again Appendix A, 
that is 

d 



^-p{x,t) = tD'-" 



hoo 



p{x,t) + / w{x — x') p{x' ,t) dx' 



oo 



(3.16) 



with p(x,0^) = 6{x). The equivalence of the two forms (3.15) and (3.16) 
is easily proved in the Fourier-Laplace domain by multiplying both sides of 
Eq. (3.14) with the factor s^"^. 

We note that the choice (i) may be considered as a limit of the choice (ii) 
as /3 = 1. In fact, in this limit we find H{s) = 1 so H{t) = t-'^/T{0) = 6{t) 
(according to a formal representation of the Dirac generalized function |13j). 
so that Eqs. (3.6)-(3.7) reduce to (3.10)-(3.11), respectively. In this case the 
order of the Caputo derivative reduces to 1 and that of the R-L derivative to 
0, whereas the Mittag-LefHer waiting time law reduces to the exponential. 

In the sequel we will formally unite the choices (i) and (ii) by defining what 
we call the Mittag-Leffler memory function 



6{t) , if /? = 1 



/f^^^(t)= rXl^' 0</5<l' (3.17) 



whose Laplace transform is 

H^^^{s) = s^-\ 0</3<l. (3.18) 

Thus we will consider the whole range < /? < 1 by extending the Mittag- 
LefHer waiting time law in (3.13) to include the exponential law (3.9). 
Remark: Equation (3.7) clearly may be supplemented by an arbitrary 
initial probability density p{x,0) = f{x). The corresponding replacement 
of 6{x)) by f{x) in (3.1) then requires in (3.4) multiplication of the 
term (1 — (f){s))/s by /(k) and in (3.6) replacement of the LHS by 
H{s) sp{k, s) — /(k) . With p{x, 0) = 6{x) we obtain in p{x, t) the 
fundamental solution of (3.7) 



4 Manipulations: rescaling and respeeding 

We now consider two types of manipulations on the CTRW by acting on its 
governing equation (3.7) in its Laplace-Fourier representation (3.6). 
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(A) : rescaling the waiting time, hence the whole time axis; 

(B) : respeeding the process. 

(A) means change of the unit of time (measurement). We replace the 
random waiting time T by a waiting time tT, with the positive rescaling 
factor T. Our idea is to take < r <C 1 in order to bring into near sight the 
distant future. In a moderate span of time we will so have a large number 
of jump events. For r > we get the rescaled waiting time density 

(f)T.{t) = (f){t/T)/T , hence 4>r{s) =4){ts) . (4.1) 

By decorating also the density p with an index r we obtain the rescaled 
integral equation of the CTRW in the Laplace-Fourier domain as 

Hr{s) sp^{K, s)-l = [w{k) - 1] p^{k, s) , (4.2) 

where, in analogy to (3.5), 

= i^itf). (4.3) 

S (p{TS) 

d 

(B) means multiplying the quantity representing -g^p{x, t) by a factor 

1/a, where a > is the respeeding factor, a > 1 means acceleration, 
< a < 1 means deceleration. In the Laplace-Fourier representation this 
means multiplying the RHS of Eq. (3.6) by the factor a since the expression 

Q 

sp{k, s) — 1 corresponds to —p{x,t). 

We now chose to consider the procedures of rescaling and respeeding in their 
combination so that the equation in the transformed domain of the rescaled 
and respeeded process has the form 

Hris) Sp^„{K, s) - 1 =a [w{k) - 1] Pr^ail^, s) , (4.4) 

Clearly, the two manipulations can be discussed separately: the choice {r > 
0, a = 1} means pure rescaling, the choice {r = 1, a > 0} means pure 
respeeding of the original process. In the special case r = 1 we only respeed 
the original system; if < r ^ 1 we can counteract the compression effected 
by rescaling to again obtain a moderate number of events in a moderate span 
of time by respeeding (decelerating) with < a <C 1. These vague notions 
will become clear as soon as we consider power law waiting times. 
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Defining now 

:=:I^ = 1Z|M. (4.5) 

as(p{TS) 

we finally get, in analogy to (3.6), the equation 

Hr,a{s) Sp^,a{l^, s) - 1 = [w{k) - 1] s) . (4.6) 

What is the combined effect of rescaling and respeeding on the waiting time 
density? 

In analogy to (3.5) and taking account of (4.5) we find 

1 + sHr,a{s) 1 - 4>{ts) 

i- ~\~ S ^ 

as ^{ts) 

and so, for the deformation of the waiting time density, the essential formula 

frAs)= ^ . (4.8) 

1 — (1 — a)(p[TS) 

Remark: The formula (4.8) has the same structure as the thinning formula 
(2.5) by identification of a with q. In both problems wc have a rescalcd 
process defined by a time scale r, and we send the relevant factors r, a 
and q to zero under a proper relationship. However in the thinning theory 
the relevant independent parameter going to is that of thinning (actually 
respeeding) whereas in the present problem it is the rescaling parameter r. 

5 Power laws and asymptotic universality of the 
Mittag-Leffler waiting time density 

We have essentially two different situations for the waiting time distribution 
according to its first moment (the expectation value) being finite or infinite. 
In other words we assume for the waiting time pdf 4>{t) either 



/•oo 

p - t' (l){t') dt' <oo, labelled as /? = 1 , 
Jo 



(5.1) 



or 

c 



(p{t) ~ ci"('^+^) for i ^ oo hence ^(i) 0</3<l, OO. 

(5.2) 
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For convenience we have dispensed in (5.2) with decorating by a slowly 
varying function at infinity the asymptotic power law. Then, by the standard 
Tauberian theory (see \12\ I68j) the above conditions (5.1)-(5.2) mean in the 
Laplace domain the (comprehensive) asymptotic form 

^(s) = 1 _ As^ + o(s^) for s ^ 0+ , < ;3 < 1 , (5.3) 

where we have 

X = p, if /3 = 1; A = cr(-/3)= / -4^,ifO</3<l. (5.4) 

T[p + 1) sm(p7r) 

Then, fixing s as required by the continuity theorem of probability theory 
for Laplace transforms, taking 

a = Ar^ , (5.5) 

and sending t to zero, we obtain in the limit the Mittag-Leffler waiting time 
law. In fact, Eqs. (4.8) and (5.3) imply as r — > with < /3 < 1, 

~ Xt" [I- Xt^s^ + oir^s^)] 1 

'Ar,Ar/3 (^) - 1 _ (1 _ _ ^ ^ 1 + ^Z? ' ^^'^f 

the Laplace transform of cf)^^^{i), see (1.1) and Appendix C. This formula 
expresses the asymptotic universality of the Mittag-Leffler waiting 
time law that includes the exponential law for /? = 1. It can easily 
be generalized to the case of power laws decorated with slowly varying 
functions, thereby using the Tauberian theory by Karamata (see again 

mm)- 

Comment: The formula (5.6) says that our general power law waiting time 
density is gradually deformed into the Mittag-Leffler waiting time density 
as r tends to zero. 

Remark: Let us stress here the distinguished character of the Mittag-Leffler 
waiting time density (f)^^^{t) = ——Ep{—t^) defined in (1.1). Considering 
its Laplace transform 

4>'''Hs) = -r^ , 0</3<l, (5.7) 
1 -|- s'^ 

we can easily prove the identity 

(^f/'(s) = ^^^^(rs/a^/^) for all r > 0, a>0. (5.8) 
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Note that Eq. (5.8) states the self- similarity of the combined operation 
rescaling-respeeding for the Mittag-Leffler waiting time density. In fact, (5.8) 
impUes 0:^'^^(t) = (f)^^{t/c)/c with c = t /a}-/^ , which means replacing the 
random waiting time T^^^ by cT^^. As a consequences, choosing a = 
we have 

0f,^,(s) = 0^'^(s) for all r>0. (5.9) 

Hence the Mittag-Leffler waiting time density is invariant against combined 
resettling with r and respeeding with a = t^. 

Observing (5.6) we can say that (j)^^{t) is a r — > attractor for any power 
law waiting time (5.2) under simultaneous rescaling with r and respeeding 
with a = Xt^. In other words, this attraction property of the Mittag- 
Leffler probability distribution with respect to power law waiting times (with 
< /3 < 1) is a kind of analogy to the attraction of sums of power law jump 
distributions by stable distributions. 



6 Passage to the diffusion limit in space 

We have again two different situations for the jump-width distribution but 
according to its second moment being finite or infinite. In other words we 
assume for the jump-width probability density wlx) (assumed for simplicity 
to be symmetric: w{x) = w{—x)) either 

/+00 
w{x) dx < OO , labelled as a = 2 , (6-1) 
'OO 

or 

w{x) r^b\x\-^°'^^^ for \x\^oo, < a < 2, 6>0. (6.2) 
Then we have the asymptotic relation, compare e.g. with |16l 120^ [2T \ I22j. 

iv{K) = 1 — fi\K\" + o{\k\") for K — > , (6-3) 

where 

/i=^ifa = 2, fi= . 77^ if < a < 2, (6.4) 

2 1 (a + 1) sm (Qvr/z) 

The above asymptotic relations are known in the framework of the attraction 
properties of the stable densities. We note that the classical book by 
Gnedenko and Kolmogorov [14J has unfortunately the wrong constant /x 
for < a < 2. As before we dispense with the possible decoration of the 
relevant power law by a slowly varying function. 
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By another respeeding, in fact an acceleration, we can pass over to space- 
time fractional diffusion processes. For this we have three choices: 

(a) : difFusion limit in space only, for general waiting time, 

(b) : difFusion limit in space only, for ML w^aiting time, 

(c) : joint limit in time and space (with power laws in both) with 
scaling relation. 

Note hat (b) is just a special case of (a) but of particular relevance (as 
we shall see). In all three cases wc rescale the jump density by a factor 
h > 0, replacing the random jumps X by hX. This means changing the 
unit of measurement in space from 1 to 1/h, with < /i ^ 1, so bringing 

into near sight the far-away space. Wc get the rcscaled jump density as 

Wh{x) = w{x/h)/h, corresponding to Wh{K) = WfiiH;). 

Choice (a): difFusion limit in space only, with a general waiting 
time law. 

Starting from the Eq. (3.6), the Laplace- Fourier representation of the 
CTRW equation, without special assumption on the waiting time density, 
we fix the Fourier variable k and accelerate the spatially rescaled process 
by the respeeding factor l/(///i"), arriving at the equation (using qh as new 
dependent variable) 

H{s) [sq^{n, s)-l\= ^ q„{K, s) . (6.5) 

Then, fixing k as required by the continuity theorem of probability theory 
for Fourier transforms, and sending h to zero we get, noting that [w{hK) — 
l]/(///i") ""I'^l") and writing u in place of qo, 

H{s) su{k,s) — 1 = u{k, s) , (6.6) 

where we still have, consistently with (3.5), 

S(j){s) (f){s) 

being (J){t) the original waiting time density. In physical space-time we have 
the integro-pseudo-differential equation 

^ H{t - t') ^u{x, t') dt' = ^D^ u{x, t), < a < 2 , (6.7) 

with — as the symbol of the Riesz pseudo-differential operator ^Dq 
usually referred to as the Riesz fractional derivative of order a, see Appendix 
B. 
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Comments: By this rescaling and acceleration the jumps become smaller 
and smaller, their number in a given span of time larger and larger, the 
waiting times between jumps smaller and smaller. In the limit there are no 
waiting times anymore, the original waiting time density (p{t) is now only 
spiritual, but still determines via H{t) the memory of the process. Eq. (6.7) 
offers a great variety of diffusion processes with memory depending on the 
choice of the function H{t). 

Choice (b): diffusion limit in space only, with a Mittag-LefHer 
waiting time law. 

We now choose in Eq. (6.7) the Mittag-Leffler memory function (3.17), 
namely 

H{t) = H^^^{t), hence -S(s) = l/s^^"^) , 
corresponding to the Mittag-Leffier waiting time law 

consistently with the time fractional Kolmogorov-Feller equation (3.15), that 
includes for /? = 1 the classical Kolmogorov-Feller equation (3.11). As a 
consequence of our spatial diffusion limit, compare with [23l I44| . we so arrive 
immediately at the space-time fractional diffusion equation 

tD^ u{x, t) = ^D^ u{x, t), 0<a<2,0</3<l. (6.8) 



Choice (c): diffusion limit combined in time and space. 

Assuming the behaviour for the waiting time density as in Eqs. (5.1)-(5.2), 
and for the jump-width density as in Eqs. (6.1)-(6.2), rescaling as described 
the waiting times and the jumps by factors r and h, starting from (4.4), 
decelerating by a factor A in time, then accelerating for space by a factor 
l/{^h°'), we obtain (compare to Section 4, case (B)), fixing s and k and 
setting, for convenience 

a{T,h) = ^, (6.9) 



Hris) 



SPr,a{T,h)[l^,- 



1 



a{T,h) [Wh{K) - 1] Pr,a{T,h){'^, 



with Wh{K) = w{hK) and 



(6.10) 



Hr{s) 



[TS] 



(rs) 



(6.11) 
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Fixing a{T,h) to the constant value 1, which means introducing the 
relationship of well-scaledness 



a{T, h) 



(6.12) 



= 1, 



between the rescaling of time and space, we get 



for T — > . 



(6.13) 



Because of ^ 

w{hn) — 1 



, for ^ ^ , 



(6.14) 



we finally get the limiting equation 



su{n,s) — \ =— |k|"u(k, s) 



(6.15) 



corresponding to Eq. (6.8), the space-time fractional diffusion equation. 

Comments on some mathematical and physical aspects 

{a) The Mittag-Leffler waiting time (choice (b)), obeying the power law 
asymptotics (5.2) with A = 1 leads from (6.7) directly to the space-time 
fractional diffusion equation (6.8), without requirement of rescaling and 
deceleration in time, and with these procedures we arrive likewise at (6.8). 
This strange fact is caused by the invariance of the Mittag-LefHer density 
to the combined effects of rescaling by r and deceleration by t^^ , expressed 
in eq. (5.9). 

Going again through our preceding deductions, we observe that the 
combined (well-scaled) passage of r and h, under the relation (6.12), towards 
zero can be split in two distinct ways into two separate passages. First wav : 
keep h fixed letting r tend to zero, then in the resulting model send also h to 
zero. Second way : interchange the order played by h and r in the first way. 
Under our power law assumptions we can transform (3.7), the basic integral 
equation of CTRW, into Eq. (3.15) (time fractional CTRW) by rescaling- 
respeeding manipulation only in the time variable, and then by rescaling 
in space followed by an acceleration into (6.8), the space-time fractional 
diffusion equation. Or we can transform (3.7) by rescaling in space followed 
by an acceleration into Eq. (6.7) (general space fractional diffusion with 
memory), and then by by rescaling-respeeding in the time variable arrive at 



(6.8). 
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(7) Where have the waiting times gone in the space-time fractional diffusion 
equation (6.8)? We can answer this question by interpreting eq. (6.10) under 
the scaHng relation (6.12) as the Laplace- Fourier representation 

= [whin) - 1] Pr,i{i^, s) , (6.16) 

of our original CTRW (3.1), whose Laplace- Fourier representation (3.6) 
coincides with (6.16) if there we delete all decorations with indices. Thus 
eq. (6.16) represents the same physical process as (3.1) but expressed in 
terms of new units 1/r and l//i of time and space, respectively. However, 
the respeeding factor a(r, h) being fixed to 1, there is no change of physical 
speed. When these new units are made smaller and smaller, moderate spans 
of time and space become numerically smaller and smaller, shrinking towards 
zero as rand h tend to zero, and likewise the waiting times and the jump 
widths shrink to zero. The distant future and the far-away space come 
numerically into near sight. As long as r and h are positive, we always 
have the same physical process, only measured in other units. The finally 
resulting space-time fractional diffusion process (6.8) remembers the power 
laws for waiting times and jumps in form of the orders (3 and a of fractional 
differentiation. 

(5) An objection could be raised against the somewhat mystical actions of 
respeeding. Namely, if the respeeding factor a in eq. (4.4) differs from 1, the 
underlying renewal process and consequently the whole CTRW are distorted. 
However, for the CTRW we carry out the actions of deceleration and 
acceleration in either order in succession or simultaneously in combination, 
and by our special choice of these factors they cancel each other in effect, 
so that there remains no physical distortion. This is particularly obvious in 
our choice (c), see the above comment (7). 

(e) Let us finally point out an advantage of splitting the passages r ^ and 
/i — s- 0. Whereas by the combined passage as in choice (c), if done in the well- 
scaled way (6.12), the mystical concept of respeeding can be avoided, there 
arises the question of correct use of the continuity theorems of probability. 
There is one continuity theorem for the Laplace transform, one for the 
Fourier transform, see [12] . Possible doubts whether their simultaneous use 
is legitimate vanish by applying them in succession, as in our two splitting 
methods. 

Discussion on the involved stochastic processes 

In our investigations we have met four types of spatially one-dimensional 
stochastic processes for the sojourn probability density p{x,t) or u{x,t). 



Hr{s) i(k,s)-1 
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For the reader's convenience let us give a list of these processes in physical 
coordinates, referring to the preceding text for details, and remind briefly 
how they can be connected by appropriate scaling and passages to the limit. 
Let us note that in all these processes the initial condition 6{x) for p(x,0+) 
or u{x,0'^) can be replaced by a more general probability density function 

m. 

(I) The integral equation for the CTRW is, see (3.1) with (3.2), 

ft r r+OD 

p{x,t) =p{x,0'^)'^{t) + (j){t-t') / w{x -x')p{x',t')dx' dt' 

Jo U-oo 

is equivalent, by the introduction of the memory function H{t), see (3.5), to 
the generalized Kolmogorov-Feller equation, see (3.7), 

/■* d /'+°° 

J H{t-t')—p{x,t')dt' = -p{x,t) + J w{x-x')p{x',t)dx' . 

(II) The time fractional Kolmogorov-Feller equation, see (3.15), 

/ + 00 
w{x - x')p{x',t)dx' , 0<P<1. 
-oo 

(III) The integro-pseudo-differential equation of space fractional diffusion 
with general memory, see (6.7), 



/' 

^0 



t Q 

Hit - t') Q^u{x, t') dt' = ^D^ u{x, t), < a < 2 . 



(IV) The space-time fractional diffusion equation, see (6.8), 

tD^u{x,t) = ^D^u{x,t) 0</?<l, 0<a<2. 



We now sketch shortly how these four evolution equations are connected in 
our theory. Eq. (I) goes over in eq. (II), likewise eq. (Ill) in eq. (IV) by the 
special choice H[t) = H^^(t) for the memory function, see (3.17). Under 
our power law assumption for the waiting time, see (5.1) and (5.2), these 
transitions can be achieved asymptotically by manipulation via rescaling 
and respeeding of the underlying renewal process. Under our power law 
assumption for the jumps, see (6.1)and (6.2), the transition from eq. (I) to 
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eq. (Ill) and from (II) to (IV) can be achieved asymptotically by passage 
to the diffusion limit only in space. Under our power law assumption for 
time and space there is a direct way from eq. (I) to eq. (IV), namely the 
well-scaled passage to the diffusion limit, for which the condition (6.12) is 
relevant. 



7 The time fractional drift process 

It is instructive to study the spatial transition to the diffusion limit for the 
Mittag-Leffler renewal process. As said in Section 3 this renewal process, 
viewed as a CTRW by treating its counting number N as a spatial variable 
X, is obtained by choosing w{x) = 6{x — 1) as the jump width density, see 
Eq. (3.3). Its waiting time density is, see (1.1), (3.13), 

<P(t) = <t>^'^{t) = -j^Ep{-t^), 0</3<l. 



We have H{s) = ^, w{k) = e"^, hence 



Sp{K, S) 



1 



(e*" - l) p(k, s) 



(7.1) 



Rescaling in space by a factor h and accelerating (because of w{Hi) = e**^ = 
1 + + o{k) for K — > 0) this pure renewal process by the factor 1/h we get 
a process 



sqh{K.,s) - 1 



1 qhiK,s), 



which as /i — > and k fixed gives 



su{k,s) — 1 =iKu{K,s), 



(7.2) 



which implies 



u{k, s) 



— IK 

We note that Eq. (7.2) corresponds to the time fractional drift equation 



(7.3) 



d 



tD^ u{x,t) = -—u{x,t) , u{x,0) = d{x), xGR+,tGR+. (7.4) 
By using the known scaling rules for the Fourier and Laplace transforms. 



f{ax) ^ a-^ f{K/a), a>0, f{bt)^b-^f{s/b), b > 
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we infer directly from (7. 3) (thus without inverting the two transforms) the 
following scaling property of the (fundamental) solution 

u{ax , bt) = h^^u{ax/h^ , t) . 

Consequently, introducing the similarity variable x/t^ , we can write 

u{x,t) = t-^U{x/t^), (7.5) 

where U{x) = u{x, 1) . 

To determine the solution in the space-time domain we can follow two 

alternative strategies related to the different order in carrying out the 
inversion of the Fourier-Laplace transforms in (7. 3). Indeed we can 

(51) : invert the Fourier transform getting u{x,s), and then invert this 
Laplace transform, 

(52) : invert the Laplace transform getting u{n, t) , and then invert this 

Fourier transform. 

Strategy (SI): Recalling the Fourier transform pair, 

ae-^^e(x) ^ 6>0, 

b — IK 

where Q{x) denotes the unit step Heaviside function, we get 

u{x,s) = sf^-^e-^^'^ &{x). (7.6) 

In view of the fact that exp(— s'^) is the Laplace transform of the extremal 
unilateral stable density of order /?, L'^^{t) (see for notation Appendix B), 
we recognize that the solution in the space-time domain can be expressed 
in terms of a fractional integral (see Appendix A) of such density, namely 
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Working in the Laplace domain we can note that the fundamental solution 
of our fractional drift equation (7.4) is simply related to that of the time 
fractional diffusion-wave equation 

tDfu{x,t) = ^u{x,t), 0<^<1, xeR, tGR+, (7.8) 

equipped with the initial conditions u{x,Q^) = 6{x) if < /3 < 1 and 
d 

—u{x,Q'^) = if 1/2 < /3 < 1. In fact, the solution of (7.8) turns out the 
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half of the solution (7.7) of our time fractional drift equation (7.4), extended 
in a symmetric way to all of R, as can be seen by factorizing eq. (7.8) as 



Indeed eq. (7.8) was solved by using the Laplace transform strategy by 
Mainardi in the 1990's, see e.g. [67\ \5E[ B5] where the reader can find 
mathematical details of the proof and instructive plots of the fundamental 
solution. Then, based on Mainardi's analysis, we can state that the required 
solution of eq. (7.4) reads 

u{x, t) = t-f^ Mpix/tl^) Q{x) , (7.9) 

where M denotes the function of Wright type defined in the complex plane 



n=0 rt=l 

(7.10) 

The M function is a special case of the Wright function defined by the series 
representation, valid in the whole complex plane, 

oo „ 

Indeed, we recognize 

M^(z) = $_^,i_^(-z) , 0</3<l. (7.12) 

Originally, Wright introduced and investigated this function with the 
restriction A > in a series of notes starting from 1933 in the framework of 
the asymptotic theory of partitions. Only later, in 1940, he considered the 
case —1 < A < 0. We note that in the handbook of the Bateman Project 
[TU] (see Vol. 3, Ch. 18), presumably for a misprint, A is restricted to be 
non negative. For further mathematical details on the M- Wright function 
we recommend [T71 [18l [23] . 

For our time fractional drift equation (7.4) we note the particular case /3 = 
1/2 for which we obtain 

P = l/2: u{x,t) = ^exp[-x^/{4t)] , x>0, t>0. (7.13) 



\n— 1 
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In the limiting case /3 = 1 we recover the rightward pure drift, 

(3 = 1: u{x,t) = ^6{x/t-l) =6{x-t), x>0, t>0. (7.14) 

In view of the fact that that M- Wright function of order /3 is related to 
the extremal unilateral stable density of order (3, see |42) . we conclude by 
displaying the alternative form of the solution of the time fractional drift 
equation: 

n(x,t) = ^^L/(t/xV/^), (7.15) 

which, compared with (7.7), shows the effect of the fractional integral on 
the stable density function L^^ . 

Strategy (S2): Recalling the Laplace transform pair, see e.g. [TOl [TO]. 

Ef.ict") ^ K(s)>|c|l/^ 
sP — c 

we get 

u{K,t) = EfsiiKt'^) , (7.16) 

from which 

1 r+oo 

u{x,t) = —VP e~^^^ E3{iKt'^)dK, (7.17) 

where VP denotes the Cauchy principal value. Because, see [TOj, Vol. 3, 
Chapter XVIII on Miscellaneous Functions, Section 18.1 Eq. (7), 

~ ^ m y ^ ±oo if < /? < 1 , 

r(i - P)y 

we see that Epliy) does not tend to zero fast enough for the integral (7.17) 
to exist as a regular improper Riemann integral. But there should be no 
problem for existence as a Cauchy principal value integral. It can be shown 
that the present strategy based on Fourier integral (7.17) provides the result 
(7.9). 

Remark: Not wanting to overload our paper we have deliberately avoided 
the concept of subordination in fractional diffusion. But, referring to [25], 
let us say that if in (7.15) we replace x by we get the subordinator, i.e. the 
probability law for generating the operational time i* from the physical time 
t, see eq. (5.20) in [25], and, in other notation, [46]. Because of its relation 
(7.16) via Fourier transform to the Mittag-Leffler function with imaginary 
argument, the probability law governing the process (7.15) sometimes is 
called the Mittag-Leffler distribution, see e.g. [46j. Although so named it 
must not be confused with our Mittag-Leffler waiting time distribution whose 
density is given by (1.1). 
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8 Conclusions 



The basic role of the Mittag-Leffler waiting time probabihty density in time 
fractional continuous time random walk (CTRW) has become well known by 
the fundamental paper of 1995 by Hilfer and Anton [29j. Earlier in the theory 
of thinning (rarefaction) of a renewal process under power law assumptions, 
see the 1968 book by Gnedenko and Kovalenko [I5], this density had been 
found as limit density by a combination of thinning followed by rescaling 
of time and imposing a proper relation between the rescaling factor and 
the thinning parameter. Likewise one arrives at this law when wanting to 
construct a certain special class of anomalous random walks, see the 1985 
paper by Balakrishnan [1], the anomaly defined by growth of the second 
moment of the sojourn probability density like a power of time with exponent 
between and 1. Balakrishnan's paper, having appeared a few years before 
the fundamental paper of 1989 by Schneider and Wyss [61], is difficult to 
read as it is written in a style different from the present one, so we will 
here not go into details. But let it be said that by well-scaled passage to 
the limit from CTRW (again under suitable power law assumptions in space 
and time) he obtained the space-time fractional diffusion equation in form 
of an equivalent integro-differential equation. Unfortunately, Balakrishnan's 
paper did not find the attention it would have deserved. However, due to the 
sad fact that the Mittag-Leffler function too long played a rather neglected 
role in treatises on special functions Balakrishnan as well as Gnedenko and 
Kovalenko contented themselves with presenting their results only in the 
Laplace transform domain; they did not identify their limit density as a 
Mittag-Leffler type function. 

Having worked ourselves for some time on questions of well-scaled passage to 
the diffusion limit from continuous time random walks to fractional diffusion, 
see [211 Ea ESI ESI ig |60] , we got from the theory of thinning the idea that it 
should be possible to carry out the passages to the limit separately in space 
and in time. In time this can be done by a combination of re-scaling time and 
respeeding the underlying renewal process (formally treating it as a CTRW 
with unit steps in space). In fact, thinning in the sense of Gnedenko and 
Kovalenko transforms the original renewal process into one that is running 
more slowly and this effect can be balanced by proper choice of the rescaling 
factor. The result of our combination of rescaling and respeeding for a 
CTRW governed by a given renewal process with a generic power law waiting 
time law is a time fractional CTRW. By another rescaling in space (now 
under power law assumption for the jumps) which can be interpreted as a 
second respeeding we arrive at the already classical space-time fractional 



24 



diffusion equation. In this way we shed new hght on the long time and wide 
space behaviour of continuous time random walks. 

In a series of comments at the end of Section 6, we have explained how, by 
what we call well-scaled passage to the diffusion limit, the transition from 
the CTRW to the space-time fractional diffusion process actually can be 
obtained by merely rescaling time and space without any respeeding at all. 
However, the separate passages to the limit are more satisfying with respect 
to mathematical rigour. 

Finally, in Section 7, we have treated the time fractional drift process as 
a properly scaled limit of the counting function of a pure renewal process 
governed by a waiting time law of Mittag-Leffler type. Our trick in finding 
the limiting waiting time law of this renewal process consists in treating it as 
a CTRW with positive jumps of size 1 so that its counting number acts as a 
spatial variable. Then, by suitably rescaling this spatial variable, we obtain 
as an interesting side result the long time behaviour of the Mittag-Leffler 
renewal process. 

Appendix A: The time fractional derivatives 

For a sufficiently well-behaved function f{t) {t > 0) we define the Caputo 
time fractional derivative of order (3 with < /? < 1 through 



C {^Df /(t); s) = / f{s) - /(0+) , /(0+) := ^lim f{t) , {A.l) 



Such operator has been referred to as the Caputo fractional derivative since 
it was introduced by Caputo in the late 1960's for modelling the energy 
dissipation in the rheology of the Earth, see [U |5] . Soon later this derivative 
was adopted by Caputo and Mainardi in the framework of the linear theory 
of viscoelasticity, see [6] . 

The reader should observe that the Caputo fractional derivative differs from 
the usual Riemann-Liouville (R-L) fractional derivative 



so that 




< /3 < 1. 




D^fit): 



d_ r 1 /(r) dT ' 

Jt [r(i-/3) io 



< /? < 1. 



(^.3) 
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Both derivatives are related to the Riemann Liouville (R-L) fractional 
integral that is defined for any order /? > as 

so that C {t f{t); s} = s'f^ f{s). Incidentally t J" t = < 7"+'^ for a,p> 
0. Then, in virtue of eqs (A.2)-(A.4), the two fractional derivatives read: 

tDf^ := tV^J^-f^ , </?<!, (A5) 

tD^ := tJ^-" tD\ 0</3<l. (A6) 

In particular, the R-L derivative of order /3 is the left inverse of the 
corresponding R-L fractional integral in that tD^ tJ^ fit) = f{t). 
We note the relationships between the two fractional derivatives (when both 
of them exist), for < /3 < 1, 

fit) = tD^ [fit) - /(0+)] = tD^ fit) - ^.^^ /(0+) . (^.7) 

As a consequence we can interpret the Caputo derivative as a sort of 
regularization of the R-L derivative as soon as /(O"*") is finite; in this 
sense such fractional derivative was independently introduced in 1968 by 
Dzherbashyan and Nersesian [9], as pointed out by Kochubei, see [SS', [34], 
In this respect the regularized fractional derivative is sometimes referred to 
as the Caputo- Dzherbashyan derivative. 

We observe the different behaviour of the two fractional derivatives (A. 2), 
(A. 3) at the end points of the parameter interval (0, 1), as it can be noted 
from their definitions in operational terms (A. 5), (A. 6). In fact, whereas for 
/? ^ 1~ both derivatives reduce to tZ)^, due to the fact that the operator 
tJ^ = I commutes with tD^ , for /3 — > O'^ we have 

^ I tD^fit) - tJ' tD' fit) = fit) - /(0+) . ^^-'^ 

The above behaviours have induced us to keep for the Riemann-Liouville 
derivative the same symbolic notation as for the standard derivative of 
integer order, while for the Caputo derivative to decorate the corresponding 
symbol with subscript *. 

For the R-L derivative the Laplace transform reads for < (3 < 1 

C {tD^ fit); s] = / 7(s) -5(0+) , 5(0+) = ^lim git) , git) := tJ^'-^^ fit) . 

iA.9) 
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Thus the rule (A. 9) is more cumbersome to be used than (A.l) since it 
requires the initial value of an extra function g{t) related to the given f(t) 
through a fractional integral. However, when /^O"*") is finite we recognize 
g(0+) = 0. 

In the limit /3 — > 1~ both derivatives reduce to the derivative of the first 
order so we recover the corresponding standard formula for the Laplace 
transform: 

£{tD'f{ty,s}=sf{s)-f{0+). (AlO) 

We conclude this Appendix noting that in a proper way both derivatives 
can be generalized for any order /3 > 1, see e.g. [T9| [53]. 



Appendix B: The space fractional derivatives 

Let us first recall that a generic linear pseudo-differential operator A, 
acting with respect to the variable x G R, is defined through its Fourier 
representation, namely 

/ + 00 
e'''''Af{x)dx = A{K)f{K), kGR (B.l) 
-oo 

where A{k) is referred to as the symbol of A, formally given as 

= (Ae"*'^^) e+^^^. {B.2) 

The fractional Riesz derivative xDq is defined as the pseudo-differential 
operator with symbol — ■ This means that for a sufficiently well-behaved 
(generalized) function /(x) (x G R) we have 

J^{.xD^f{x);K} = -\Krf{K), KGR. (S.3) 

The symbol of the Riesz fractional derivative is nothing but the logarithm 
of the characteristic function of the generic symmetric stable (in the Levy 
sense) probability density, see (HJHSIE^. Noting —\k\" = — (k^)"/^, we 
recognize that 

'"^ = -[-1^) ■ f^-^' 

In other words, the Riesz derivative is a symmetric fractional generalization 
of the second derivative to orders less than 2. In an explicit way the Riesz 
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derivative reads, for < a < 2, 
^D" f(x) = S^f(x) 

= r(i+a)^^y^ e^- 

where in the L.H.S we have also adopted the alternative and illuminating 
notation introduced by Zaslavsky, see e.g. [56]. This operator is referred 
to as the Riesz fractional derivative since it is obtained from the inversion 
of the fractional integral originally introduced by Marcel Riesz in the late 
1940's, known as the Riesz potential, see e.g. [57]. It is based on a suitable 
regularization of a hyper-singular integral, according to a method formerly 
introduced by Marchaud in 1927. 

Remark: Straightforward generalization to the Riesz-Feller derivative of 
order a and skewness 9 is possible. Such pseudo-differential operator is 
denoted by us as 

^D^ , with < a < 2 , 61 e R, \e\ < min{a, 2 - a} . (-B.6) 
In this case we have 

r{,D1f{x)-K] = -\K\^i^'''^''^f{K) KGR. {B.l) 
In an explicit way the Riesz-Feller derivative reads, for < a < 2, 



.D^ m - ^ {sin Ka + .)./2] f ^^^^^ 



(i?.8) 



Notice in (B.7) that i^siguK ^ exp[i (sign k) 6* 7r/2]. Thus the symbol of 
the Riesz-Feller fractional derivative is the logarithm of the characteristic 
function of the more general (strictly) stable probability density, closely 
following the Feller parameterization, see [11^ 112] revisited by the present 
authors in [20]. According to our notation, the strictly stable density of 
order a and skewness 6 is denoted by L^{x). We note that the allowed 
region for the parameters a and 9 turns out to be a diamond in the plane 
{a, 6} with vertices in the points (0, 0), (1, 1), (2, 0), (1, —1), that we call the 
Feller- Takayasu diamond, see Fig. 1. For more details we refer the reader to 
[42] , where series representations and numerical plots of the stable densities 
L^(x) are found. In particular, we recall that the extremal stable densities 
obtained for 9 = ita with < a < 1 are unilateral, with support in R"^, 
respectively. 
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Fig. 1 The Feller- Takayasu diamond 



Appendix C: The Mittag-LefHer function 

The Mittag-LefHer function with parameter /3 is defined as 

oo „ 

It is an entire function of order (3 and reduces for /? = 1 to exp(2;) . For 
detailed information on the functions of Mittag-LefHer type the reader may 
consult e.g. [101 [El ESI [39l ESj [57] and references therein. 
Hereafter, we find it convenient to summarize the most relevant features of 
the functions 

^{t):=Ef,{-t^), 0</3<l, {C.2) 

^{t) := -jEpi-t") , </?<!, {C.2,) 

that turn out to be the most relevant functions of Mittag-LefHer type for 
our purposes. Both of them reduce to the exponential function exp(— t) in 
the limit as /3 — > 1. 

We begin to quote their expansions in power series of (convergent for 
t > 0) and their asymptotic representations for t ^ oo, 

iJ r(/3n + l) vr ' ^ ' 



n=0 



n=0 ^ ' 
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The Laplace transforms of ^{t) and (p{t) can easily be obtained by 
transforming the series (C.4), (C.5) term by term, respectively: they read 

$(.) = ^, ^{s) = -^, Rs>0. (C.6) 

For < (3 < 1 both functions ^(t), (p{t) keep the complete monotonicity 
of the limiting exponential function of exp(— t) . Complete monotonicity of 

a function f{t) means, for n = 0,1,2,..., and t > 0, ( — 1)" — ^ /(O — 
or equivalently, its representability as (real) Laplace transform of a non- 
negative function or measure, see e.g. [T2] . 

Recalling the theory of the Mittag-Leffler functions of order less than 1, we 
obtain for < /? < 1 the following representations, see e.g. |19j . 



^r(t) = / — — — dr, t>0, (C.7) 

7T Jo r2/3 + 2r/3 cos(/37r) + 1 ' - ' ^ ^ 

, sin(/37r) r^Q-rt 

Mt) = / — -— dr, t>0. (C.8) 

' vr Vo r2/3 + 2r/3 cos(/37r) + 1 ^ ' 

In Figs 2 and 3 we exhibit plots of the functions ^{t) and (j){t), respectively 
in logarithmic and linear scales. 




Fig. 2 The functions ^{t) (left) and (j){t) (right) in logarithmic scales 




Fig. 3 The functions ^{t) (left) and (j){t) (right) in linear scales 
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